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, Abstract. For many classically chaotic systems it is believed that 

the quantum wave functions become uniformly distributed, that is 
the matrix elements of smooth observables tend to the phase space 
average of the observable. In this paper we study the fluctuations 
of the matrix elements for the desymmetrized quantum cat map. 
f-H ' We present a conjecture for the distribution of the normalized ma- 

, trix elements, namely that their distribution is that of a certain 

weighted sum of traces of independent matrices in SU{2). This 
is in contrast to generic chaotic systems where the distribution 
^2 ' is expected to be Gaussian. We compute the second and fourth 

moment of the normalized matrix elements and obtain agreement 
with our conjecture. 

m 
> 

(N 

Q I A fundamental feature of quantum wave functions of classically chaotic 

■ systems is that the matrix elements of smooth observables tend to the 

phase space average of the observable, at least in the sense of conver- 
'-^ [ gence in the mean 121 ^] or in the mean square jIH]. In many 

J3 systems it is believed that in fact all matrix elements converge to the 

micro-canonical average, however this has only been demonstrated for 
a couple of arithmetic systems: For "quantum cat maps" jTU], and 
I conditional on the Generalized Riemann Hypothesis^ also for the mod- 

I ular domain [12], in both cases assuming that the systems are desym- 

metrized by taking into account the action of "Hecke operators" . 

As for the approach to the limit, it is expected that the fluctuations 
of the matrix elements about their limit are Gaussian with variance 
given by classical correlations of the observable [71IS]. In this note we 
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study these fluctuations for tlie quantum cat map. Our finding is that 
for this system, the picture is very different. 

We recall the basic setup [HI 13 IH HHI (see section |21 for further back- 
ground and any unexplained notation): The classical mechanical sys- 
tem is the iteration of a linear hyperbolic map A G SL{2, Z) of the 
torus = R^/Z^ (a "cat map"). The quantum system is given by 
specifying an integer N, which plays the role of the inverse Planck con- 
stant. In what follows, will be restricted to be a prime. The space 
of quantum states of the system is Hn = L'^{Z/NZ). Let / G C°°It^) 
be a smooth, real valued observable and Op^(/) : Hn Hn its quan- 
tization. The quantization of the classical map A is a unitary map 
Un{A) of Hn. 

In [TO] we introduced Hecke operators, a group of commuting uni- 
tary maps of Hn, which commute with Un{A). The space Hn has an 
orthonormal basis consisting of joint eigenvectors {'ipj}jLi of Un{A), 
which we call Hecke eigenf unctions. The matrix elements (Op^(/)'?/'j, ipj) 
converge^ to the phase-space average J^a f{x)dx [TU]. Our goal is to 
understand their fluctuations around their limiting value. 

Our main result is to present a conjecture for the limiting distribution 
of the normalized matrix elements 

F^^ ■■= ViV ((Op^(/)^„^,) - f{x)dx 

For this purpose, define a binary quadratic form associated to A by 

Q{x, y) = cx^ + ((i — a)xy — hy^ , A = 

For an observable / G C°°(T^) and an integer u, set 

n=(ni,n2)eZ2 
Q(n)=u 

where f{n) are the Fourier coefficients of /. 

Conjecture 1. As N ^ oo through primes, the limiting distribution 
of the normalized matrix elements Fj^'^ is that of the random variable 

where are independently chosen random matrices in SU (2) endowed 
with Haar probability measure. 

^For arbitrary eigenfunctions, that is ones which are not Hecke eigenfunctions, 
this need not hold, see 
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This conjecture predicts a radical departure from the Gaussian fluc- 
tuations expected to hold for generic systems [7| E] . Our first result 
confirms this conjecture for the variance of these normalized matrix 
elements. 

Theorem 2. ^5 ^ oo through primes, the variance of the normal- 
ized matrix elements Fj^^ is given by 

(1.1) ^j:\^r\'-^m)=j:\f*(^)\'- 

j=l u^O 

For a comparison with the variance expected for the case of generic 
systems, see Section IHTTl A similar departure from this behaviour of the 
variance was observed recently by Luo and Sarnak ^12] for the modular 
domain. For another analogy with that case, see section 

We also compute the fourth moment of F^^^ and find agreement 
with Conjecture^ 

Theorem 3. The fourth moment of the normalized matrix elements is 
given by 

j=l uj^O 

as N —>■ 00 through primes. 

Acknowledgements: We thank Peter Sarnak for discussions on his 
work with Wenzhi Luo [T^ . 



2. Background 

The full details on the cat map and its quantization can be found 
in pni- For the reader's convenience we briefly recall the setup: The 
classical dynamics are given by a hyperbolic linear map A G SL{2, Z) 
so that a; = (g) G T^ I— > Ax is a symplectic map of the torus. Given 
an observable / G C°°(T^), the classical evolution defined by A is 
f f o A, where (/ o A){x) = f{Ax). 

For doing quantum mechanics on the torus, one takes Planck's con- 
stant to be and as the Hilbert space of states one takes TYat : = 
L^(Z/A^Z), where the inner product is given by 

Q mod N 
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The basic observables are given by the operators Ti^{n), n G Z^, 
acting on^p e L^{Z/NZ) via: 

(2.1) (T^K, n,)ij) (Q) = e^e(^)^(Q + n,). 

where e(x) = e^'^*^'. 

For any smooth classical observable / G C°°(T^) with Fourier ex- 
pansion /(x) = J2nez^ /(^)e(nx), its quantization is given by 

Op;v(/) := E n^)TN{n) . 

2.1. Quantum Dynamics: For A which satisfies a certain parity con- 
dition, we can assign unitary operators Uj\f{A), acting on L^(Z/A^Z), 
having the following important properties: 

• "Exact Egorov": For all observables / G C°°(T2) 

UNiA)-' OpM)UNiA) = Op^(/ o A). 

• The quantization depends only on A modulo 2N: If A = B 
mod 2N then Un{A) = Un{B). 

• The quantization is multiplicative: if A, B are congruent to the 
identity matrix modulo 4 (resp., 2) if is even (resp., odd), 
then [ini ESI 

Un{AB) = Un{A)Un{B) 

2.2. Hecke eigenfunctions. Let a, be the eigenvalues of A. Since 
A is hyperbolic, a is a unit in the real quadratic field K = Q(a). Let 
D = Z[«], which is an order of K. Let v = (fi,f2) G be a vector 

such that vA = av. If A = ^ , we may take v = {c, a — a) . Let 

/ := Z[fi, ^2] = Z[c, a — a] C D. Then / is an D-ideal, and the matrix 
of a acting on / by multiplication in the basis vi,V2 is precisely A. The 
choice of basis of / gives an identification 1 = 7? and the action of D 
on the ideal / by multiplication gives a ring homomorphism 

Mat2(Z) 

with the property that det(t(/5)) = A/'(/9), where M : Q(a) — > Q is the 
norm map. 

Let C{2N) be the elements of 0/2N0 with norm congruent to 1 
mod 2N, and which congruent to 1 modulo 40 (resp., 20) if N is even 
(resp., odd). Reducing t modulo 2N gives a map 

i2N ■■ C{2N) SL2{Z/2NZ). 
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Since C{2N) is commutative, the multiplicativity of our quantization 
implies that 

{UN{i2Nm : /3 e C} 
forms a family of commuting operators. Analogously with modular 
forms, we call these Hecke operators, and functions ip ^T-Cn that are si- 
multaneous eigenfunctions of all the Hecke operators are denoted Hecke 
eigenf unctions. Note that a Hecke eigenfunction is an eigenfunction of 
UN{i2N{a)) = Un{A). 

The matrix elements are invariant under the Hecke operators: 

(Op^(/)^„^i) = (Op^(/o5)^„^,), B e C{2N) 

This follows from ipj being eigenfunctions of the Hecke operators C{2N). 
In particular, taking /(x) = e{nx) we see that 

(2.2) {TNin)^P,,ij,) = {TN{nB)^j,^j) . 

2.3. The quadratic form associated to A: We define a binary qua- 
dratic form associated to A = 

Q{x, y) = cx^ + (d — a)xy — hy^ 

This, up to sign, is the quadratic form A/'(a;c+?/(a— a))/A/'(/) induced 
by the norm form on the ideal J = Z[c, a — a] described in Section 
where M{T) = ^D/I. Indeed, since I = Z[c,a — a] and O = Z[l,a] 
we have A/'(/) = |c|. A computation shows that the norm form is then 
sign(c)Q(x,y). 

By virtue of the definition of Q as a norm form, we see that A 
and the Hecke operators are isometrics of Q, and since they have unit 
norm they actually land in the special orthogonal group of Q. That is 
we find that under the above identifications, C{2N) is identified with 
{B e SO{Q, Z/2NZ) ■.B = I mod 2}. 

2.4. A rewriting of the matrix elements. We now show that when 
ip is a Hecke eigenfunction, the matrix elements {Op]sf{f)ip,ip) have a 
modified Fourier series expansion which incorporates some extra in- 
variance properties. 

Lemma 4. If m,n G are such that Q{m) = Q{n), then for all 
sufficiently large primes N we have m = uB mod for some B G 
SO{Q,Z/NZ). 

Proof. We may clearly assume Q{m) ^ because otherwise m = n = 
since Q is anisotropic over the rationals. We take A^ a sufficiently 
large odd prime so that Q is non-degenerate over the field Z/NZ. If 
A^ > then Q{m) ^ mod A^ and then the assertion reduces 
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to the fact that if Q is a non-degenerate binary quadratic form over 
the finite field Z/A^Z (A^ 7^ 2 prime) then the special orthogonal group 
SO{Q,Z/NZ) acts transitively on the hyperbolas {Q{n) = z/}, z/ 7^ 
mod N. □ 

Lemma 5. Fix m,n & T? such that Q{m) = Q{n). If N is a suffi- 
ciently large odd prime and ip a Hecke eigenf unction, then 

Proof. For ease of notation, set e(n) := (— l)"^"'^. By Lemma |3] it 
suffices to show that if m = uB mod for some B e SO{Q, Z/NZ) 
then e{n){TN{n)ip,ip) = e{m){TN{m)tp,tp) . 
By the Chinese Remainder Theorem, 

SO{Q, Z/2NZ) ~ SO{Q, Z/NZ) x SO{Q, Z/2Z) 

(recall N is odd) and so 

C{2N) ~ {5 G S0{QZ/2NZ) : B = I mod 2} ~ SO{Q, Z/NZ) x { J} 

Thus if m = nB mod for B G SO{Q, Z/NZ) then there is a unique 
B G C(2A^) so that m = nB mod A^. 

We note that e{n)TN{n) has period A^, rather than merely 2A^ for 
TnIu) as would follow from ()2.H) . Then since m = nB mod A^, 

e{m)T]s[{m) = e{nB)TN{nB) = e{n)T]s[{nB) 

(recall that B G C{2N) preserves parity: nB = n mod 2, so e{nB) = 
€{n)). Thus for ip a Hecke eigenfunction, 

t{m){TN{m)i),i)) = e{n){TN{nB)i),i:) = t{n){TN{n)i),i)) 

the last equality by □ 

Define for z/ G Z 

f*{v) := 

neZ2;Q(n)=iy 

and 

(2.3) K(V^) := yiV(-l)"^"MT;v(ri)V',V') 

where n G Z^ is a vector with Q{n) = v (if it exists) and set (■?/') = 
otherwise. By LemmaElthis is well-defined, that is independent of the 
choice of n. Then we have 

Proposition 6. Iftpisa Hecke eigenfunction, f a trigonometric poly- 
nomial, and N > A"o(/), then 
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To simplify the arguments, in what follows we will restrict ourself to 
dealing with observables that are trigonometric polynomials. 



3. Ergodic averaging 

We relate mixed moments of matrix coefficients to traces of certain 
averages of the observables: Let 

1 



(3.1) 



D n) 



\Ci2N)\ 



BeC{2N) 



The following shows that D{n) is essentially diagonal when expressed 
in the Hecke eigenbasis. 

Lemma 7. Let D be the matrix obtained when expressing D(n) in 
terms of the Hecke eigenbasis {'4'i}f=i- If N is inert in K, then D is 
diagonal. If N splits in K, then D has the form 
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D21 


D22 
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D = 
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V 
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• Dnn/ 



where ipi,ip2 correspond to the quadratic character of C{2N). More- 
over, in the split case, we have 



A 



for 1 < i, j < 2. 



Proof. If N is inert, then the Weil representation is multiplicity free 
when restricted to C{2N) (see Lemma 4 in [5].) If is split, then 
C(2N) is isomorphic to [Ta/NTa)* and the trivial character occurs with 
multiplicity one, the quadratic character occurs with multiplicity two, 
and all other characters occur with multiplicity one (see [n], section 
4.1). This explains the shape of D. 

As for the bound on in the split case, it suffices to take f{x,y) = 
^^nix^ny^ for some ni,n2 G Z. We may give an explicit construction of 
the Hecke eigenfunctions as follows (see jJJ, section 4 for more details): 
there exists M G SL2{Zi/2NZ) such that the eigenfunctions ipi,ip2 can 
be written as 



i^i = VN-UN{M)So, i)2 



N 



N 



f/^(M)(l-5o) 
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where 5o(x) = 1 if a; = mod A^, and 6o{x) = otherwise. Setting 
01 = V^^o and 02 = wzji'^ ~ ^o), exact Egorov gives 

Dij = {TNiini,n2))i>i,i>j) = (TAr((n'i, n2))0i, 0^) 

where (^'i,?t-2) = {ni,n2)M mod A^. Since we may assume n not to 
be an eigenvector of A modulo A^, we have n[ ^ mod and n'2 ^ 
mod A^. Hence 

= {TN{{n[,n',))<P^,<P,) = e{^)6o{0 + n[) = 

since n'l ^ mod A^. The other estimates are analogous. □ 

Remark: In the split case, it is still true that D^j <^ N~^^'^ for all 
but this requires the Riemann hypothesis for curves, whereas the 
above is elementary. 

Lemma 8. Let {ipi}fLi be a Hecke basis o/Hn, and let k,l,m,n G Z^. 
Then 

N 

Y,{TN{m)ij„ V;.)(T^(n)^„^,) = tr {D{m)D*{n)) + 0{N-') 
1=1 

Moreover, 

N 
i=l 

= tr {D{k)D*{l)D{m)D*{n)) + 0{N-^) 
Proof. By definition 

N N 



^{TNim)'4ji,ipi){TNin)ipi,ip,i) = ^ D{m)iiD{n)ii 
1=1 1=1 

On the other hand, by lemma Q 

AT 

tr {D{m)D{ny) = Duim)D^) + D2i{m)D^) + ^ Du{m)Du{n) 

i=l 

where Di2{m), D2i{m), Di2{n) and D2i{n) are all O^N^^^"^). Thus 

TV 

J2{TN{m)iJi, ilJi){TN{n)tlj,,tlJi) = tr {D{m)D{ny) + 0{N-^) 

i=l 

The proof of the second assertion is similar. □ 
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4. Proof of Theorem El 

In order to prove Theorem El it suffices, by Proposition IHl to show 
that as ^ cxD, 

^E^^(^.)m) - E(trf/.trf/,) = 1^1^ _^ 

where f/^, U,^ G SU2 are random matrices in SU2, independent if u ^ fi. 

Proposition 9. Let {4'i}fLi be a Hecke basis o/Hn- If N > NQ{ji,u) 
is prime and fi,u ^ mod A^, then 

iV ) otherwise. 

Proof. Choose m^n such that Q{m) = /i and Q{n) = v. By fj2.3p 
and Lemma IHl we find that 

^ N N 

= (_l)'™i'"2+"i"2tr {D{n)D{my) +0{N-^) 
By definition of D{n) we have 

D(n)D(m)* = ^ 5^ T;v(n5i)T^(mi?2)* • 

We now take the trace of both sides and apply the following easily 
checked identity (see (EH)), vahd for odd N and 5i, ^2 G C{2N): 



tT{TN{nBi)TN{mB2y) = 
We get 

1 ^ 

^ -^^mim2+nin2 



-l)™i"^2+nin2jv if nEi = mod N, 

otherwise. 



\C(2N)\^ 

nB\=mB2 mod TV 

^ \{B eC{2N) ■.n = mB mod N}\ + 0{N-^] 



\C{2N)\ 

which, since |C(2A^)| = ± 1, equals 1 + 0(iV"^) if there exists B e 
C{2N) such that n = mi? mod A^, and 0{N^^) otherwise. Finally, for 
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sufficiently large (i.e., > Nolfi, z/)), Lemma ^ gives that n = rriB 
mod for some B G C{2N) is equivalent to n = u. □ 

5. Proof of theorem El 
5.1. Reduction. In order to prove Theorem |21 it suffices to show that 



(5.1) ^^K(^,)Va(^,)V;(^,)K(V',) ^E(tr^7,tr[/Atr[/^tr[/,) 
i=i 

where U^, Ux, and f/^ are independent random matrices in SU2. 

Let S* C Z"^ be the set of four-tuples (k, A, fi, v) such that k = A, /i = 
u, or K = fi, X = u, or K, = u, X = fi, but not n = X = fi = u. 

Proposition 10. Let be a Hecke basis ofTiN OL^d let k, X, fi,^ & 

Z. If N is a sufficiently large prime, then 

N (2 + 0{N-') zfK = X = fx = u, 

■'=1 [0(A^"^/^) otherwise. 

Given Proposition it is straightforward to deduce ()5.H) . we need 
only to note that E((tr?7)^) = 2, E((trf/)2) = 1, and E(trf/) = 0. 

The proof of PropositionlTUlwill occupy the remainder of this section. 
For the reader's convenience, here is a brief outline: 

(1) Express the left hand side of ()5.1|) an exponential sum. 

(2) Show that the exponential sum is quite small unless pairwise 
equality of k. A, /x, z/ occurs, in which case the exponential sum 
is given by the number of solutions (modulo A^) of a certain 
equation. 

(3) Determine the number of solutions. 
5.2. Ergodic averaging. 

Lemma 11. Choose k, l,m,n G Z^ such that Q{k) = n, Q{1) = A, Q{m) = 
yU, and Q{n) = v. Then 



\C{2N)\* 



^ n A^ 

Bi,B2,S3,B4GC{Af) 
kBi-lB2+mB3-nB4,=0 mod N 



t{uj{kBi, -IB2) + cj(m53, -^54)) 
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The proof of Lemma ITD is an extension of the arguments proving the 
analogous (j4.1j) in the proof of Proposition El and is left to the reader. 

5.3. Exponential sums over curves. In order to show that there is 
quite a bit of cancellation in ()5.2j) when pairwise equality of norms do 
not hold, we will need some results on exponential sums over curves. 
Let X be a projective curve of degree di defined over the finite field 
Fp, embedded in n-dimensional projective space P" over Fp. Further, 
let R{Xi, . . . , Xn+i) be a homogeneous rational function in P", defined 
over Fp, and let d2 be the degree of its numerator. Define 



where a is the trace from Fpm to Fp, and the accent in the summation 
means that the poles of R{x) are excluded. 

Theorem 12 (Bombieri yy. Theorem 6). If did2 < p and R is not 

constant on any component T of X then 



In order to apply Bombieri's Theorem we need to show that the 
components of a certain algebraic set are at most one dimensional, and 
in order to do this we show that the number of points defined over Fat 
is 0{N). (Such a bound can not hold for all if there are components 
of dimension two or higher.) 

Lemma 13. Let a,b ^ ¥^[0]. If a 7^ and the equation 

7i = «72 + b, 7i, 72 e C{N) 
is satisfied for more than two values 0/72, then 6 = and N'{a) = 1. 

Proof. Taking norms, we obtain 1 = A/'(a) +A/'(6) +tr(a672) and hence 
tr(a672) is constant. If ab 7^ 0, this means that the coordinates {x,y) 
of 72, when regarding 72 as an element of F^, lies on some line. On 
the other hand, N{pi2) = 1 corresponds to 72 satisfying some quadratic 
equation, hence the intersection can be at most two points. (In fact, 
we may identify C{N) with the solutions to — Dy"^ = 1 for x, y G Fat, 
and some fixed D G F^f.) □ 

Lemma 14. Fix k,l,m,n G and let X be the set of solutions to 



k-lB2 + niB-i -nB4 = mod A^, B2, B3, B4 G C{N) 

If Q{k),Q{l),Q{m),Q{n) ^ mod A^, then \X\ < 3(A^ + 1) for N 
sufficiently large. 




SUR,X)\ < {dl + W2 - 3rfiK/' + dl 
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Proof. We use the identification of tlie action of C{N) on witli tlie 
action of C{N) on F7v[a]. The equation 

k-lB2 + niB-i - nBi = mod N 

is then equivalent to 

K - A/52 + /^/^s - l^l^A = 
where G C{N) and k, A, yU, z/ G FAr[a;]. We may rewrite this as 

and letting P' = we obtain 

If — fiP' = then n — XP2 = 0, and since Q{1), Q{m) ^ mod 
implies that A, are nonzero^, we find that (32 and (3' are uniquely 
determined, whereas /34 can be chosen arbitrarily. Thus there are at 
most |C(A^)| solutions for which v — fij3' = 0. 

Let us now bound the number of solutions when i> — fi(3' 7^ 0: after 
writing 

K, - A/?2 = /S^iiy - n(3') 

as 

+ 77,P2 = 



Lemma IT^ gives (note that k 7^ since Q(/c) ^ mod iV) that there 
can be at most two possible values of P21 Pa for each and hence there 
are at most 2|C(A^)| solutions for which u — 7^ 0. Thus, in total, 
X can have at most |C(A^)| + 2|C(A^)| < 3(A^ + 1) solutions. □ 

5.4. Counting solutions. We now determine the components of X 
on which e ^ *('^(^^i'~'^2)+^(™-^3,-".g4)) j jg constant. 

Lemma 15. Assume that Q{k), Q{1), Q{itl), Q{n) ^ mod A^, and let 
Sol(fc, /, m, n) he the number of solutions to the equations 

(5.3) kBi - IB2 + mB^ - uB^ = mod 

(5.4) uj{kBi, -IB2) + uj{mB^, -uB^) = -C mod 

where Bi G C{N). If C = mod A^ and N is sufficiently large, then 
(5.5) 

r 2\C{N)\' tfQik) = Q{1) = Q{m) = Q{n), 

Sol(fc,/,m,n) = <^ \C{N)\^ + 0{\C{N)\) if {Q{k),Q{l),Q{m),Q{n)) G S, 
lo(|C(A^)|) otherwise. 



"'Recall that Q, up to a scalar multiple, is given by the norm. 
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On the other hand, if C ^ mod then 

Sol{k,l,m,n) = 0{\C{N)\). 

Proof. For simplicity^, we will assume that is inert. It will be con- 
venient to use the language of algebraic number theory; we identify 
{Z/NZf with the finite field F^2 = ¥n{VD) by letting m = {x,y) 
correspond to /i = x + yy/D. First we note that if n = (z, w) corre- 
sponds to u then 

u!{m, n) = xw — zy = Im((a; + yVD){z + wVD)) 

where Im(a + b\/D) = b, and hence uj{m,n) = Im(/7z/). 

Thus, with {k,l,m,n) corresponding to (z/i, 2^2, t's, 1^4), the values of 
Q{k) , Q{1) , Q{m) , Q{n) modulo are (up to a scalar multiple) given 
by Ar(z/i),Ar(z/2),Ar(z/3),Ar(z/4). Putting /i, = for A e C{N), we 
find that uj{kBi, —IB2) + uj{mB3, —nB^) = —C can be written as 

Im(/Zr/i2 + Ji3fJ"i) = c. 

Now, kBi — IB2 + mi?3 — uBi = mod is equivalent to yUi — /i2 = 
f^A — fJ'3- Taking norms, we obtain 

A/'(/ii) + Mifi2) - tr(7ir/i2) = Mifii) + MifJ^s) - trifMl^s) 

and hence 

tiijMfis) = tr(7ir/i2) + N^ + N3-Ni-N2 
if we let A^'j = N{ui). Since tr(/i) = 2Re(/i) = 2Re(7l), we find that 

2 Re(7Ii/i4) = 2 Re(;Ui7I^) + Ni + N^-Ni-N2 
On the other hand, Im(7ii/i2 + J^f^A) = C implies that 

Im(7I^/i4) = - Im(7J7/i2) + C = lm{jjijl^) + C 

and thus 

7^3/^4 = /il7i2 + K 

where i^' = (A^4 + A^s - A^i - A^2)/2 + C\/D. Hence we can rewrite flO) 
and ()5.4p as 

{Af3A*4 = /^iAi2 + K 
/^i + yUs = + /i4 
/i, = z/i/3,, AeC(Ar)fori = l,2,3,4. 

^The split case is similar except for possibility of zero divisors, but these do not 
occur when k, I, m, n are fixed and A'' is large enough. 
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Case 1 {K ^ 0). Since /ij = z/j/5j with Pi G C{N), we can rewrite 

7Ii/i4 = fiijt^ + K 



as 



and hence 



f34/f33 = ^{^iV^f3i/l32 + K). 
Z/3Z/4 



Applying lemma IT^ with 71 = P^/Ps and 72 = /9i//92 gives that P1/P2, 
and hence /ii/i2, must take one of two values, say Ci or C2. But /ii/i2 = 
Ci implies that /xi = /^2^ and hence = /X3 '"^'^ • We thus obtain 

/^2(1 - ^) = /ii - /^2 = - /is = /^sll 

Now, if /xi 7^ /i2 then both 1 — and 1 — '"^'^ are nonzero. Thus 
/i2 is determined by ^3, which in turn gives that /ii as well as /i4 are 
determined by /is. Hence, there can be at most C{N) solutions for 
which /xi 7^ /i2- (The case fiijl2 = C2 is handled in the same way.) 
On the other hand, for /ii = /i2 we have the family of solutions 

(5.6) /ii = /i2, /i4 = /i3 

(note that this implies that C = Im(/i7/i2 + /i3/i4) = 0.) 
Case 2 {K = 0). Since K = and /ii = yU2 + Ai4 — /^s we have 
7Zi/i4 = yUi/^ + K = {n2 + H4- /^s) A*2 

and hence 

/i4(7¥ - 7^) = (/U2 - Ats)/^ 

If /i2 — /^s = 0, we must have /xi = /i4, and we obtain the family of 
solutions 

(5.7) /i2 = /^s, Ail = /i4 

On the other hand, if /i2 — /^s 7^ 0, we can express /i4 in terms of /i2 
and yUs: 

_ Ai2 - Ai3 _ A^2 - /i2/W3 

/i4 — = =Ai2 — 77 = /^S, 

Ai3 - At2 ^3 - /i2/i3 

which in turn gives that 

(5.8) /^l = /i2 + /i4 - /i3 = /i2 + = =/U2-A'3 

/is - /i2 

/i2 - /is . s , /i2 - /is /i2 - /is /i27i3 - iVs 

= = — =(/is - /i2) + = — =/i2 = = — =/is = —= rr/i2 

/is - /i2 /is - /i2 /^S - /i2 /i2/iS - A'2 
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Summary. If K there can be at most 2\C{N)\ "spurious" solu- 
tions for which /ii ^ ^2] other than that, we must have 

On the other hand, if X = 0, then either 

or 

1^2 - AtS ^2 - /^/^3 Ai2 - Ai3 At2Ai3 " ^3 

IJ'3 -1^2 Ns- /X2A*3 A*3 - A*2 A*2A*3 " ^2 

We note that the first case can only happen if A'"! = N2 and A^3 = iV4, 
the second only if N2 — and Ni — N4, and the third only if N2 — N4 
and Ni = Ns- Moreover, in all three cases, C = Im(ii') = Im(7Zr/^2 + 
Jisfj^^) = 0. We also note that if N2 = N3, then the third case simplifies 
to III = 112 and /X3 = We thus obtain the following: 

If C 7^ then X 7^ and there can be at most 0{N) "spurious 
solutions" . 

If C = and Ni = N2 = N^^ N4 then X = and the solutions are 
given by the two families 

/^2 = /^3, /^i = /^4 

and 

N2 - Jl^Hs fi2'p3 - N3 

If C = and Ni ^ N4 ^ N2 ^ then X = and there is a family 
of solutions given by 

/^2 = 1^3, 1^1 = 1^4- 

Similarly, if C = and Ni = N3 ^ N2 ^ then = and there 
is a family of solutions given by 

/"2 - fki 1^2-/^3 

At4 = = =A*2, fJ'l = = =A*3 

/^3 — f^2 f^3 — 1^2 

If C = and Ni ^ N2 ^ ^ N4 then X 7^ 0, in which case we 
have a family of solutions given by 

III = II2, IJ'3 = 1^4 

as well as 0{N) "spurious" solutions. 

Finally, if C = and pairwisc equality of norms do not hold, then 
we must have K {if K = then 7J3/X4 = 1111I2 + K implies that 
N3N4 = N1N2, which together with Ni + N2 = N3 + N4 gives that 
either AT-^ = iVg, iVs = N4 or Ni = N4,N2 = N^) and in this case there 
can be at most 0{N) "spurious" solutions. 
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Now Lemma 13] gives that pairwise equality of norms modulo N im- 
plies pairwise equality of Q{k), Q{1), Q{m), Q{n). □ 

5.5. Conclusion. We may now evaluate the exponential sum in (j5.2j) 

Proposition 16. If Q{k),Q{l),Q{m),Q{n) ^ mod then, for N 
sufficiently large, we have 

(5.9) 

't{uj{kBi, -IB2) + cj(mS3, -nBi)) 



E 



e 



Bi,B2,B-i,Bi^C(N) 
kBi-lB2+m.B3-nB4=0 mod TV 



N 



r 2|C(iV)|2 + 0{\C{N)\) tfQik) = Qil) = Q{m) = Q{n), 
= I |C(iV)|2 + 0(|C(iV)|) if (Qik), Qil), Qim), Q(n)) e S, 
[0{\C{N)f/^) otherwise. 

Proof. Since both uj{kBi, —IB2) + uj{mB^, —uB^ and kBi — IB2 + 
mi?3 — ni?4 are invariant under the substitution [Bi, B2, B^, B4) —>■ 
(B'Bi, B'B2, B'B^, B'Bi) for B' e C{N), we may rewrite the left hand 
side of (El) as |C(A^)| times 



(5.10) ^{ 

B2,B3,Bi&C{N) ^ 
k-lB2+mB3-nBi=0 mod N 



t{uj{k, -IB2) + i^{mB^, -uBi)) 
N 



Let X be the set of solutions to 

k-lB2 + niB-i -nBi = mod A^, B2, B3, B^ e C{N). 

By Lemma El the dimension of any irreducible component of X is 
at most 1. The contribution from the zero dimensional components 
of X is at most 0(|C(A^)|). As for the one dimensional components, 
Lemma IT^ gives that uj{k, —IB2) + uj{mB^, —uB^ cannot be constant 
on any component unless pairwise equality of norms holds. Thus, if 
pairwise equality of norms does not hold, Bombieri's Theorem gives 
that (ICTH) is 0{N^'^) = 0(|C(A^)|i/2). 

On the other hand, if uj{kBi, —IB2) + uj{mB-i, —uB^) equals some 
constant C modulo A^ on some one dimensional component, then LemmalT^ 
gives the following: C = mod A^, and ()5.10p equals Sol(A;, /, m, ra), 
which in turn equals |C(A^)p or 2|C(A^)p depending on whether Q(A;) = 
Q{1) = Q{m) = Q{n) mod A^ or not. □ 

Proposition now follows from Lemma ^2 and Proposition on 
recalling that \C{N)\ = \C{2N)\ =N±1. 
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6. Discussion 



6.1. Comparison with generic systems. It is interesting to com- 
pare our result for the variance with the predicted answer for generic 
systems (see 13 EJ), which is 



(6.1) 



oo „ 

V / fo{x)MAhc)dx 



where fo = f — /rp2 f{y)dy- Using the Fourier expansion and collecting 
together frequencies n lying in the same A-orbit this equals 



E 



f=-oo 05>^ngZ2 mg{Z2-0)/(A) nGm(A) 

where {A) denotes the group generated by A. We can further rewrite 
this expression into a form closer to our formula (jl.lj) by noticing that 
the expression e{n) := (— l)'^i"2 is an invariant of the A-orbit: e(n) = 
e{nA), because we assume that A = I mod 2. Thus we can write the 
generic variance (|6.1|) as 

2 



(6.2) 



E 

mg(Z2-0)/(A> 



E (-l)"^"V(n) 

n(im{A) 



m 



The comparison with with our answer ^^^q '^Q{n)=ui~^)^^"^^ fi^) 
(ll.lj) . is now clear: Both expressions would coincide if each hyperbola 
G : Q{n) = u} consisted of a single A-orbit. It is true that each 
hyperbola consists of a finite number of A-orbits for u 0, but that 
number varies with u. 

6.2. A differential operator. There is yet another analogy with the 
modular domain, pointed out to us by Peter Sarnak: We define a 
differential operator L on C°°(T^) by 

1^ = -— Q(-,-) 

47r^ dp' dq 

so that Lf{n) = Q{n)f{n). 

Given observables f,g, we define a bilinear form B{f,g) by 

B{f,g) = J2f*i^)9*ii^) 

SO that (cf. Conjecture HI) B{f,g) = E(XjXg) and by Theorem 
B{f, f) is the variance of the normalized matrix elements. 
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It is easy to check that L is self adjoint with respect to B, i.e., 
B{Lf,g) = B{f,Lg). Note that L is also self-adjoint with respect 
to the bilinear form derived from the expected variance for generic 
systems fl6.1|) . ()6.2j) . This feature was first observed for the modular 
domain, where the role of L is played by the Casimir operator [12^ (c.f. 
Appendix 5 of Sarnak's survey [MJ). 

6.3. Connection with character sums. Conjecture ^ is related to 
the value distributions of certain character sums, at least in the case 
of split primes, that is primes for which the cat map A is diagonal- 
izable modulo N. Let M e SL2{Z/2NZ) be such that A = MDM-^ 
mod 2N. In |TT] we explained that in that case, all but one of the nor- 
malized Hecke eigenfunctions are given in terms of the Dirichlet char- 
acters X modulo as tp^^ := ^J^^Un{M)x- We can then write the 

matrix elements {TN{n)i/j^, ip^) as characters sums: Setting {1711,1712) = 
nM, we have 

Q mod N 

and Conjecture [T] gives a prediction for the value distribution of these 
sums as x varies. 
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